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Abstract 

Making use of the extended Ginzburg Landau theory, which includes the 
fourth order derivative term, we study the vortex state in a magnetic field 
parallel to the c axis. The vortex core structure is distorted due to this higher 
order term, which reveals the fourfold symmetry. Further, this distortion 
gives rise to the core interaction energy, which favors a square lattice tilted 
by 45° from the a axis. The critical field of this transition is determined. The 
magnetization diverges at the transition. This suggests the transition is of 
the first order. 
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I. INTRODUCTION 



After a few years of controversy, d-wave superconductivity in the hole-doped high 
cuprates appears to be finally established However, the electron-doped high Tc cuprates 
appear to be described by s-wave superconductivity 0,^. 

D-wave superconductivity manifests itself as fourfold symmetry of the vortex state when 
a magnetic field is applied either parallel to the c axis or within the a-b plane [^]. In particular 
to the study of the vortex lattice in the vicinity of the upper critical field |^ and the quasi- 
particle spectrum around a single vortex [[7|,^ in a magnetic field parallel to the c axis indicate 
that the square vortex lattice tilted by 45° from the a axis should be most stable except 
in the immediate vicinity of the superconducting transition temperature Tc. Indeed such a 
square lattice, though elongate in the a direction has been seen in YBCO monocrystals by 
small angle neutron scattering (SANS) and scanning tunneling microscopy (STM) 
at low temperature and in a low magnetic field. On the other hand, the fourfold symmetry 
predicted for the density of states near the vortex core appeared not to have been seen 
by STM [0] in YBCO monocrystals. This, we believe, indicates the failure of the quasi- 



classical approximation used in these theoretical analysis. Indeed, recent studies [0,|T3I of 
the Bogoliubov-de Gennes (BdG) equation clearly indicate not only the breakdown of the 
quasi-classical approximation for YBCO, but also the presence of the extended states with 
small energies (say \E\ < O.IA) which exhibits clearly the fourfold symmetry anticipated 
from the square vortex lattice. 

More recently a very similar square vortex lattice have been seen in ErNi2B2C, YNi2B2C 
and LuNi2B2C by SANS P,|l|] and in YuNi2B2C by STM imaging [|l|]. Although super- 
conductivity in borocarbides is believed to be conventional s-wave [|l^], the above square 
lattice together with the presence of antiferromagnetic phase in closely related borocarbides 
suggest strongly that superconductivity in borocarbides will be of (i-wave as well [|1^]. In- 
cidentally the square vortex lattice and related vortex lattice transition have recently been 



studied using the generalized London equation |T3,0,[T8|. The phenomenological free energy 
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used by these authors resembles the one for d-wave superconductivity. 

The object of this paper is twofold, (i) Making use of the extended Ginzburg-Landau 
equation, we first study a single vortex line in a magnetic field parallel to the c axis. Unlike 



Refs. ||Ty| , P!8[] we believe that the modification of the vortex core structure is of prime im- 
portance. Indeed, the vortex exhibits the fourfold symmetry, which will have a number of 
consequences. For example it will modify the quasi-particle spectrum around a vortex. One 
more significant fact is that this will generate vortex core interaction energy, which favors 
the alignment of two vortices either parallel to (1, 1, 0) or (1, —1,0). Indeed a similar vortex 
solution has been found numerically previously by Enomoto et. al. ||19|. But our analytical 



result is of prime importance in the following, (ii) From a study of the two-vortex problem, 
we consider the vortex lattice for a class of isoceles. We find in the low field limit {i.e. 
B ~ Hci{t)) the vortices form a triangular lattice as in a conventional s-wave superconduc- 
tor. When the magnetic field increases, the triangular lattice transforms first gradually and 
then suddenly to the square lattice when B = Her- In the temperature range not very far 
from T, {i.e. ^T^ < T < we predict 

H,r = 0.52A{-\nt)-^/'^K-^H,2{t), (1.1) 

where t = T/Tc and k is the Ginzburg-Landau parameter. Though the 5-dependence of the 
apex angle 6 we obtained is rather similar to the ones obtained in Refs. [T^JT^, the detail 
is quite different. For example we find the vortex lattice transition is of the first order in a 
sharp contrast to Ref. OJl^, where it is found of the second order. Furthe r the present 



model describes 6 dependence on B more consistent with SANS result [0 than that of Ref. 

18| , which may suggest that the core interaction between two vortices will be more critical 
than the term arising from the anisotropy of the magnetic interaction considered in Refs. 

1^,|T^. Unfortunately the related SANS study for high cuprates is not available at this 



time of writing. 
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II. EXTENDED GINZBURG-LANDAU EQUATION AND SINGLE VORTEX 

PROBLEM 



We consider a weak-coupling model for d-wave superconductors 



Extending the 



procedure used by Ren et.al. pT|, we obtain 

+2{dl-d^.?]]A(r 



16(47rT)4 
21C(3) 



5(9^ + a 



2\2 

y' 



A(r)|^A(r), 



(2.1) 



which is converted into the dimensionless form 



(l + {dl + dl) + 6 [b{dl + dlf + 2{dl - dl 



AW = lAMPAfr) 



(2.2) 



where we have introduced 



A(r) 



2_ (47rT)2(-lnt) 



2(47rT)2(-lnt)' ' ' 21C(3) 

t = T/Tci and rescaled r — ^(T)r, A(r) — * A(T)A(r). Here and dy are gauge invariant 
differential operators and we define the small parameter e = 31C(5)(— lnt)/196C(3)^ ~ 
0.114(-lnt). 

Equation ( ^.1[ ) is written down basically in |T^, though we ignore a few terms of the order 
of (— Int)^ since they are of secondary importance in what follows. Here we concentrate on 
the effect of the e-term, which is the basic symmetry breaking term. 

Assume that A(r) is given by 



A(r) = g{r)e''^ + e (e^*^a(r) + e-^''^(3{r) + 7(r)) e*^. 



(2.3) 



Substituting this in Eq. ( p.2|) we find g{r) for r ^ 1; 



q(r] = 1 r r r 

^ 2 8 16 



(2.4) 



and equations for a{r),f3{r) and 7(r) for r ^ 1; 
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where 



A{t) (dl + - ^)) a(r) = g{Tf{2a{T) + /3(r)), 

B{t) + (l + {dl + - (5{t) = ^(r)2(a(r) + 2/?(r)), 
C(r) + (l + + - ^)) 7(r) = ^(r)^37(r), 
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(2.6) 
(2.7) 



(2.8) 
(2.9) 
(2.10) 



Then we find 
5 
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-6627 - 184 c 253, ^ a 



and 
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7(rj = — 9r —r — ? 



(2.11) 
(2.12) 



(2.13) 



In this solution we find a free parameter c, which fortunately does not show up in the core 
interaction term which we are going to discuss in the following section. Note that the choice 
c = 5/4 makes the first few terms symmetric; a{r) — 5/2r~^ + (l — 11 logr)5/4r~^ • • • , /3(r) = 
— 5/2r~^ + (l + ll logr)5/4r~^ • • •. We will also discuss in the next paragraph that the choice 
c ~ 5/4 is necessary to have approximate solutions. 

For later purposes, it is convenient to introduce the interpolation expressions which give 
the correct asymptotics for r — > 0. We find 



g(r) = tanh [ 1 — cisech — ) tanh^ 1 — C2sech 

Co 2r^ V Co/ Co 8r^ V co 



r \ ^ , 9 r 
tanh — 
Co 



a{r) 



— 2 1 7 ^ f 1 \ —4 1 11 ^ 

-r tanh h t log r r tanh — 

2 C3 V4 4 ^ ; C3 

5 2 1 5 ^ f ^ 1 \ 4 1 Q ^ 

— r tanh h — I logr r tanh — 

2 C3 V4 4 ^ y C3 



(2.14) 

(2.15) 
(2.16) 



where cq = 1.71, ci = 0.80, C2 = 1.35. The way to fix these constants is the following. Using 
the GL equation ( |2.2| ), we can express all the constants Ci,C2--- by Ci. The constant ci 
can be obtained by performing numerical integration of the GL equation with the boundary 
conditions g{0) = 0, lim^^oo fi'('") = 1- In principle, we can apply the same procedure to a(r) 
and (3{r). However, we simply start from the ansatz ( p^.lSD and ( |2.16| ) which are given from 
( p.liP and (|2.12|) by introducing suitable powers of tanhr/cs, and observe that these with 
C3 = 2.5 and c = 5/4 agree very nicely with the numerical results obtained by Enomoto 
et. al. |jl9|. We show in Fig.l a{r) and j3{r) as function of r. These are compared with 
Sf[^\r) and 8/1^ ^lijr) in Enomoto et. al.. We see our analytic expressions are very close to 
the numerical ones in Enomoto et. al. We have not shown 7(r) as this term is somewhat 
different from the one in Enomoto et. al. since our starting equation is different. 

III. INTERACTION BETWEEN TWO VORTICES 

Before studying the regular vortex lattice, let us consider the two-vortex problem. We 
assume that two vortices are placed at (0, 0) and {dcosO, dsm.6) and 1). The free 

energy in dimensionless units is given by 

n = j dh (-|A|2 + |9,A|2 + |5j,A|2 

-emdl + dl) + 2{dl - dl))^\' + \\^\' + ^^') (3-1) 

where h = b{r) is the local magnetic field. Making use of the usual approximation A(r) = 
Ani/(r-ri), where 

/(r) = {g{r) + e(a(r)e^^^ + /3(r)e-^^<^ + 7(r))) e^^ (3.2) 

is the single vortex solution, g{r) ~ tanhr and neglecting 7(r) which is irrelevant for the 
fourfold symmetry, we obtain 

^ two-vortex - ~\ j (tanhr + ecos40(a(r) + /9(r)))^ 
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x(tanhr' + ecos40'(Q;(^') + 
~ (A - 2ai - 2aie{a{d) + (3{d)) cosiO) , (3.3) 

where A is the area and 

ai = y" rfV (2sechV - sechV) = y (^ln2 + ~ 6.854. (3.4) 

On the other hand the magnetic interaction between two vortices is given by |jA'o(^) 
(the London formula) where Kq{z) is the modified Bessel function. Strictly speaking the 
magnetic interaction is also modified due to the higher order term (see, for example, Ref. 
18). Indeed the correction term decays like d""^ with d, but this term does not contain the 
extra «;-dependence. Therefore the correction term to the magnetic interaction is completely 
neghgible when k ^ 1 as in high Tc cu prates. Therefore the core interaction give a strongly 
directional energy ~ (Z~^cos4^, while the magnetic energy is isotropic as in conventional 
s-wave superconductor. 



IV. VORTEX LATTICE 

Let us consider a vortex lattice where lattice points are given by r;„i — 
r;,,„(cos^;,OT,sin^/_OT) = /(i(cos6',sin6') + m(i(cos6', — sin 6*), where /,m are integers d = 
\J sin{2e)B -> ^0 quantum. For later convenience, we separate the lattice into 

even and odd lattices as v'f^ = r;^^^(cos 6'|-^^)j, sin ^^l"^) = 2/ci(cos^,sin^) + 2md(cos^, — sin^) 
and r["l = rl'^{cos 6^'^, sin e^^l) = {21 + l)(i(cos^, sin ^) + (2m + l)d(cos^, - sin^). Note / 
and m run over all possible integers. Then the free energy of the vortex lattice is given by 

n^A |A-aie'n^-el0aie'n^$:'4^cos4O +^^<^e'E%(^)> (4-1) 
^ V '."^ ^^'"^ / ^"^ V A / 

where = B/cpo the vortex density per unit area. Here we consider only the vortex core 
interaction between two vortices, since the three vortex interaction is exponentially small 
when d/^ ^ 1. Further, we have neglected the fourfold symmetric term in the magnetic 
interaction term since it is proportional to e/k^. So except for the condensation energy 



{—^A), the second term and the last term are proportional to B, while the core interaction 
energy (the third term) is proportional to B^. As the magnetic field increases from B = 
Hci{t), the third term becomes more and more dominant and for B > Her the square vortex 
lattice will be established. The last term in Eq. (|4.1|) contains the sum 



p—e.o 



l,m 

T 



'Ko ((/V' + m'^^'Y') (((/ - 1/2) + (m - l/2f^^'Y' 



Lm 



Lm 



where /i = 2d sin 6/ X, fi' = 2d cos 6 / X. Following the argument by Fetter et.al. namely, 
using the integral representation of the function Kq{x) and two Poisson summation formulas 
(see Appendix), we can rewrite these infinite summations. Then the last term in Eq. ( |4.1|) 
becomes (for A rf) 



271 ^r, ^ — f '^l,m 

T 



K 

2t\ 



Lm 



Att 1 , uu' 1 ,^ , 
h - In— (1-7) 



(4.2) 



1 t-^ , I „ I /,2 A'' I 2 A'' 



m 



Lm 



/i 



-l)'+'" + exp [-n{Pi^ + m^f] 



The angle ^min which minimizes the free energy is obtained by studying the function 



sin2 2^cos4^, 



Lm 



\H*it)J ((/ + m)2 sin^ ^ + (/ - m)2 cos2 ^) 

Lm \ 



where 



H*{t) 



98C(3)2(27r) 



1/2 



7r(/2cot 6 + m^tanl 



:^~ 5.64667(-lnt)-i/2:^, 



,155aiC(5)(-lnt)/ k 
Then the minimization of f{6) gives Fig. 2 where the apex angle Omm is shown as function 
of B/Hcr where 



Her = 0.52A{-\nty^^^K-^Hc2it). 



(4.3) 



For B > Her the square lattice is fully established. Note also that dO/dB diverges at 5 = Her 
indicating the possible phase transition. Earlier a similar 9-B curve was obtained within 
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the generalized London equation |T^,|T8|. However, the present result appears to be more 



consistent with the observed B dependence of 9 by SANS from ErNi2B2C at T = 3.5K [jl2 
Inserting 6 determined thus into Eq. ([4.1|), we find the free energy 



where the first term 



..n„ + ?^7f«, (4.4) 



cr 



no = + -^B + -^B 

2 200 /t'0o 



2vrA2 1 00 1 ■ 
^^+2^°S2^-2('-^^ 



(4.5) 



depends on in a non-singular way, and the second term is 



f^^n. f # ) ) . (4.6) 



We show f{B/Hcr) as a function of B/Hcr for < B/Hcr < 1-2 in Fig. 3. A cusp 
is observed at 5 = H^r- The magnetization —M = dQ/dB has a singularity at 5 = 
Her due to the cusp in f{B/Hcr)- Fig. 4 shows the singular part of the magnetization 
-Msingi = dJ{B/Her)/d{B/Her) for < B/H^r < 1-2. Fig.3 and Fig.4 show clearly this 
phase transition is of the first order with negative latent heat {i.e. less entropy in the square 
lattice), which should be readily accessible experimentally. Suppose we choose 9 — 90° as an 
order parameter, which does not have discontinuity. The approach to zero, however, is very 
sharp and it can not be described by power law as the usual second order phase transitions. 
The transition can be classified as a special kind of the first order transi tion, even though 
the order parameter is continuous. A detailed study on the nature of this phase transition 
is under consideration 



V. CONCLUDING REMARKS 

By analyzing the extended Ginzburg-Landau equation for rf-wave superconductor, we 
discover that the vortex core contains a long range fourfold term which is proportional to 
r"'' cos 40 when r > ^. The effect of this term on the quasi-particle spectrum is under 
current study. This fourfold term gives rise to the vortex core interaction, which favors the 

9 



orientation of two vortices parallel to the diagonal directions (1, 1, 0) and (1, —1, 0). In the 
low field regime we find that the vortex lattice transforms from triangular to square as B 
increases and that the last transition to the square lattices is very steep and of the first 
order. The present result app ears to describe very well the vortex transition observed in 
ErNi2B2C, though the superconductivity in borocarbides is believed to be s-wave. Turning 
to high Tc cuprates there is no similar measurement available even for YBCO monocrystals. 
On the other hand, if we put k, = 100, Hc2{0) = 120 Tesla for YBCO, we estimate Her = 1 
Tesla, which is consistent with the observation of the square lattice at low temperature and 
in a magnetic field of a few Tesla. Clearly a parallel measurement of the i?-dependence of 
the apex angle 6 in high cuprates is highly desirable. 

Coming back to the vortex lattice transformation in the vicinity of -B ^ Hc2{t), it is 
shown that the transition is again continuous in contrast to an earlier analysis [^]. In 
particular the full transition to the square lattice is completed at t = 0.81. Therefore it is 
now possible to draw a vortex lattice phase diagram in the T-B plane. 

We expect also that the directional core potential not only modifies the equilibrium 
vortex lattice configuration but also the collective mode, the elastic and dynamic response 
of the vortex lattice. At this moment we can say only that d-wave superconductivity should 
bring a profound change in our understanding of the vortex motion. 
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APPENDIX 



In this Appendix, we list some useful formulas for studying the free energy f2 (O) of 
the vortex lattice with the apex angle 9. 
We have to treat the lattice sums 

1,'m 

= E^o (((/ - 1/2) V + {m- l/2f^^"f 

l,m 

The Poisson sum formulas 



Eexp (-/VV4r) = ^exp (-47rV/2;x^ 

I ^^ I 



Eexp (-(/ - l/2)V/4r) = ^^("1)^^^? (-47rMV 



can be obtained from Jacobi's imaginary transformations for the elliptic theta functions; 
i^3{v,t) = e"*/V-i/2g-^™V-^3(y/^^_l/^) and i^4{v,t) = e'^/V-^/^g-Tr^V-^^lt^/r, -1/r). 
Using the argument by Fetter et.al. p2l, we obtain 



1 ( ( u u' \ expf-7r(/2ii^ + m2 4 

2. _ 



2 n\ ' 
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FIGURES 



Fig. 1 Plots of a{r) and /3(r). 

Fig. 2 Appex angle 2^min as a function of B/Hcr where 2^min — 90° and 120° correspond 
to the square lattice and the triangular lattice with hexagonal symmetry, respectively. 
Fig. 3 Singular part of the free energy / as a function of B/Hcr- 
Fig. 4 Singular part of the magnetization —Mgingi as a function of B/Hcr- 
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